HST multi-epoch imaging of the PSR B0540-69 system unveils a 
highly dynamic synchrotron nebula. 1 
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ABSTRACT 

PSR B0540-69 is the Crab twin in the Large Magellanic Cloud. Age, energetic 
and overall behaviour of the two pulsars are very similar. The same is true for 
the general appearance of their pulsar wind nebulae (PWNe). Analysis of Hubble 
Space Telescope images spanning 10 years unveiled significant variability in the 
PWN surrounding PSR B0540-69, with a hot spot moving at ~ 0.04c. Such 
behaviour, reminiscent of the variability observed in the Crab nebula along the 
counter-jet direction, may suggest an alternative scenario for the geometry of the 
system. The same data were used to assess the pulsar proper motion. The null 
displacement recorded over 10 y allowed us to set a 3a upper limit of 290 km s~ x 
to the pulsar velocity. 



1 Istituto Universitario di Studi Superiori di Pavia, Via Luino 4, 27100 Pavia, Italy 
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Subject headings: stars: neutron — pulsars: individual (PSR B0540-69) 



1. Introduction 



PSR B0540-69 in the Large Magellanic Cloud (LMC) is one of the youngest pulsars 
known to dat^l (characteristic age r ~ 1600 years) and one of the very few first observed 
at wavelengths othe r than radio. It was discovered in X-rays by the Einstein Observatory 
( Seward et al. 1984 ). soon detected to pulsate in the optical ((Middleditch Sz Pennypacker 
1985|), but detected in radio only 10 years later iManchester et al.lll993h . PSR B0540-69 
is a fast (~ 50 ms), classical pulsar with a rotational energy loss very similar to that of 
the Crab (~ 1.5 x 10 38 e rg s _1 ). Thus, the detection of a polarized plerion-like structure 



( jChanan fc Helfandlll990l ) came as no surprise, making it the most Crab-like of the Crab-like 
remnants. 



With the first high resolution optical images of the field (jCaraveo et al.l Il992l ) it was 
possible to identify the pulsar optical counterpart (V ~ 22.6), disentangling it from the 
surrounding structured plerio n. This picture was confirmed by a snapshot Hubble Space 
Telescope (HST) observation (jCaraveo et al.lll998l ) which clearly resolved, within ~ 4" from 
the pulsar, the plerion, elongated in the northeast-southwest direction. 



Using early Chandra data. iGotthelf fc Wand (120001 ) and lKaaret et al.l (120011 ) performed 
a morphological study of the plerion in the X-ray band and found a noticeable similarity 
(accounting for di f ferent distance) with the Crab pulsar- wind nebula (PWN). Furthermore, 
Gotthelf fc Wang! (120001 ) unveiled the presence of a brighter PWN region south-west of the 
pulsar. Somehow in analogy with the Crab case, they suggested such region to belong 
to a torus around the source, since it appears perpendicular to a much fainter structure 
protruding from the pulsar and tentatively identified as a jet. 



Caraveo et al.l (120001 ) performed a multiwavelength analysis of the PSR B0540-69 PWN 



morphology by superimposing Chandra and HST images, finding a good correlation between 
the optical and X-ray structures with the PWN emission enhanced in bot h cases SW of the 
pulsar, i.e. along the putative torus proposed by IGotthelf fc Wangj (120001 ) . 



More recently, detailed studies of the system, based on both narrow and wide-band HST 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc. under contract No. NAS 5-26555. 
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http: / / www.atnf.csiro.au/research/pulsar / psrcat / 
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observations (jSerafimovich et al.ll2004j ; iMorse et al.ll2006l ). further strenghtened the similarity 
with the Crab owing to the presence of a cage of filam entary ejecta possibly originating, at 
least in part, in a pre-supernova mass ejection phase (jCaraveo et al.lll998l ). The complex 
interaction between the PWN and such envelope is dominating the plerion multiwavelength 
morphology, as confirmed by iPetre et al.l (120071 ) , based on the analysis of a deep Chandra 
observations. 



Interestingly, ISerafimovich et al.l (12004 ). using two HST images taken ~ 4 years apart, 
reported a tentative pulsar proper motion measurement of 4.9 ±2.3 mas y _1 (for a pulsar 
projected velocity of 1190 ± 560 km s _1 ), aligned with the putative southern jet of the 
PWN. This would make P SR B0540 -69 the third pulsar, after the Crab ( Caraveo fe Mignani 
19991 : iNg fe Romanfeood ) and Vela JCaraveo et al.lboOll : bodson et alil2003h pulsars, with a 
proper motion aligned with its PWN jet and, possibly, with the spin axis, which would have 
important consequences for pulsar kick models as well for the studies of the pulsar/PWN 
interactions. 

In this Letter, we report the results of our analysis of the extended PSR B0540-69 
HST archived dataset which allowed us to study possible morphological changes in the PSR 
B0540-69 PWN and to assess the pulsar proper motion over a longer time baseline. 



2. HST data analysis and results 



Recent HST observations of PSR B0540-69 were performed on 2005, November 15 th 
with the Wide Field and Planetary Camera 2 (WFPC2), using the wide band F555W (480 
s) and the medium band F547M (1040 s) filters. Such data add up to WFPC2 ob servations 



collected in 1995 and in 1999 with several different filters (see IMorse et al.ll2006l . for a log 
of the observations). In order to study the PWN variability as well as the pulsar proper 
motion an accurate superposition of multi-epoch data is required. We decided as a first step 
to use only the observations performed using the F555W filter, having a better signal to 
noise. Moreover, the significant WFPC2 geometric distortion, crucial fo r a correct proper 
motio n measurement, has been accurately mapped for the F555W filter ([Anderson fe King 
20031 ). The selected dataset, which also provides a time baseline of > 10 year s, includes the 
recen t 2005 observation, as well as the original 1995, October 19 th (600 s) one (jCaraveo et al. 



1998|). 



The data were retrieved from the ST-ECF archive and reprocessed on-the-fly using 
the most appropriate reference files. Data reduction and analysis was performed using the 
IRAF/STSDAS, Midas and FTOOLS packages. Individual frames collected during each visit 
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were combined to remove cosmic rays hits, and averaged. Residual cosmic ray traces were 
removed using specific algorithms within Midas. 

In order to register the frames, we followed t he procedure we already applied in several 
previous astrometri c works with HST 



Caraveo et al. 1996; De Luca et al. 2000 



see e.g. 

Mignani et al.ll2000l ; ICaraveo et al.ll200ll ). A relative reference frame was defined for each 
image, selecting a sample of good (well resolved, not saturated, not extended, not too close 
to the CCD border) reference sources. In the crowded field of PSR B0540-69, 85 good 
sources were identified. Their position was evaluated fitting a 2-D gaussian function to their 
intensity profile, with a resulting uncertainty of 0.02-0.06 pixel per coordinate. The position 
of the pulsar optical counterpart was evaluated in the same way, with an uncertainty of order 
0.03-0.04 pixel per coordinate. The coordinates of the reference stars and of the pulsar were 
then c orrected for the WFPC2 geometric dist ortion using the mapp ing of I Anderson fc King 
J200J, as well as for the "34 th row" defect JAnderson fc KineJll999h . 



Next, the 1995 reference grid was assumed as a reference and was aligned along Right 
Ascension and Declination according to the telescope roll angle. Then, we computed the 
best plate transformation (accounting for independent shift and scale factor for each axis, 
as well as for a rotation angle) between the two grids of reference stars. We applied an 
iterative clipping routine, to discard reference stars yielding larger residuals. After rejecting 
15 stars, we obtained a very good frame superposition, with rms uncertainties of 0.06 pixel 
per coordinate. 



2.1. Proper motion of the pulsar 

The 2005 pulsar position was translated to the 1995 reference frame, to compute the 
pulsar displacement over the 10.1 year time span. Such a displacement turned out to be of 
0.01±0.09 pixel in R.A. and of 0.02±0.08 pixel in Dec, i.e. statistically null. Using the well 
calibrated WFPC2 plate scale, we set a 3a upper limit to the source proper motion of ~ 1.7 
mas yr _1 . At the known pulsar distance (51 kpc), such a limit corresponds to an upper limit 
to the total pulsar velocity (projected on the plane of the sky) of ~ 410 km s -1 . 

As a further check, we included in our astrometric analysis the WFPC2 observations 
performed through the F547M filter on 2005, November 15 th and on 1999, October 17 th 
(800 s). Indeed, the F547M filter (pivot wavelength 5483A, AA = 483A FWHM) has a 
narrower bandpass than the F555W one (pivot wavelength 5439 A, A A = 1228A FWHM), 
but the pivot wavelength is essentialy the same. Thus, the wavelength dependence of the 
geometric distortion should not induce any bias when using the correction optimized for 
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the F555W filter. The analysis was performed as above, using the same reference stars. 
The superposition accuracy to the 1995 reference grid turned out to be accurate within 
~ 0.08 pixels per coordinate, (i.e. only slightly less accurate than for the superposition of 
the F555W data), with no evidence for systematic effects. Also in this case, no significant 
displacement was measured for the pulsar. Combining F555W and F547M data yields a 
tighter 3a proper motion upper limit of 1.2 mas y _1 , corresponding to a projected velocity 
of ~ 290 km s" 1 . Such a limit, c ompu ted using a well-tested, robust algorithm, supersedes 
the result by ISerafimovich et al.l (120041 ). based on a much shorter time baseline. 



2.2. Variability of the PWN 

The epoch-to-epoch coordinate transformation was then used to rebin and superimpose 
the images, in order to search for possible variations of the PWN morphology. The resulting 
images were not corrected for the geometric distortion. However, we note that at the PWN 
position, imaged at t he center of the PC chip , the maximum distorsion correction is ~ 0.05 
pixel per coordinate (jAnderson &: King 20031 ) . i.e. small enough for our goals. 



As a first step, we checked consistency of photometry between the observations per- 
formed through the same filter at different epochs. To this aim, we compared the count 
rates from 85 reference stars using simple aperture photometry. The 2005 F555W obser- 
vations yield systematically lower count rates, with an average 2005-to-1995 ratio of 0.86 
(0.07 rms), while the 2005-to-1999 ratio for the F547M observations is 0.90 (0.08 rms). Such 
values were used to re-normalize the 2005 images. 

We started using the 1995 and 2005 F555W images, characterized by a better signal-to- 
noise. A striking change in the brightest portion of the PWN is immediately apparent when 



comparing the two images. As noticed by ICaraveo et al.l (120001 ) a definite surface brightness 



maximum is seen in the PWN, South- West to the pulsar optical counterpart. Such a feature, 
which we will call the "hot spot", is resolved (7.2 pixel FWHM, or 0.081 pc at the LMC 
distance) and lies ~ l"l away from the pulsar (0.27 pc) in the 1995 image. Inspection of the 
2005 image shows that the hot spot is displaced by ~ Of! 5 in the SW direction with respect 
to its location in 1995. This is apparent from Fig. [I], where the two images are compared. 
Fitting a gaussian function to the hot spot profile, we estimated that the feature moved 
by 0M6 ± 0'.'02, or 0.11 pc at the LMC distance, corresponding to a velocity of ~ 0.037c 
assuming simple linear motion. 

The apparent displacement of the hot spot corresponds to a very large local variation 
in the PWN surface brightness. Using a 8 x 8 pixel aperture (~ 0V36 x 0'.'36) centered on the 
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1995 hot spot position, the count rate is seen to decrease by 25 ±2% between 1995 and 2005 
(the uncertainty does not account for systematic errors in the image renormalization). We 
have computed the ratio of the 2005 to the 1995 surface brightness using the same 8x8 pixel 
aperture, selecting 200 positions to cover the whole PWN (within ~ 2'.'5 from the pulsar). 
Such an exercise proved that the region of the hot spot is by far the most active part of 
the PWN. The observed r.m.s. variability on the above 200 PWN regions is ~ 8%, which 
reduces to ~ 4% when considering only the brightest 100 regions. Of course, the possible 
systematics involved in the 2005 to 1995 renormalization do not affect such a conclusion. 

The hot spot is clearly seen in the image collected in 2005 with the F547M filter, at 
a position consistent with the one apparent in the F555W filter. The ratio of the 2005 
F547M and F555W images does not show any significant feature at the hot spot position, 
which suggests the hot spot emission to be dominated by continuum. Indeed, considering 
a 4 pixel radius aperture centered at the hot spot, we estimated the ratio of the observed 
background-subtracted count rates to be 0.4 ±0.1. This is fully consistent with an expected 
value of 0.45, evaluated with the WFPC2 ETc|, assuming a power-law spectrum of spectral 
index a —1. 



Thus, we can use the 1999 F547M image (see Fig.l of ISerafimovich et al.l 120041 ) to 
constrain the position of the hot spot at a third epoch. Results are shown in Fig. [21 The 
hot spot peak in 1999 lies ~ 0'.'38 West of its 1995 position, while in 2005 it is seen ~ 0'.'33 
South of its 1999 position. The hot spot morphology is also seen to vary, the feature being 
more extended in 1999. Such results argue against a simple outward motion of the feature 
and prove a dramatic variability of the PWN in the SW region. 

The detection of large time variability in the PWN of PSR B0540-69 makes its similarity 
with the Crab Nebula even more compelling. Thus, it seems natural to compare in some de- 
tail the optical phenomenology of the two systems. The hot spot in the PSR B0540-69 PWN 
is definitely larger and m ore distant from the puls ar than the bright, highly variable "wisps" 
seen in the Crab Nebula (IHester et al.lll995l . |2002| ). It is somewhat reminiscent (as for physi- 
cal dimensions, distance to the pulsar and temporal behaviour) o f a large, roughly ar c-shaped 
structure in the outer Crab nebula, first noticed in the optical by lHester et al.l (119951 ) because 
of its outstanding variability on a time scale of 6 years (see their Figure 12d ). Such a feature 



is also prominent and highly variable at radio wavelengths (see Fig. 2 of iBietenholz et al. 



! http: / /www. stsci.edu/hst /wfpc2 / software/ 



iSerafimovich et alJ (j200j) in their spatially- resolved study of the PWN spectrum used a region ( "area 2" ) 
encompassing the hot spot. Emission from such a region, largely contributed by the hot spot, is consistent 
with a power law of spectral index 1.6 ± 0.4. 
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2004 ). While the nature of such feature is not understood, it is almost certainly re lated to 
energy outflows in the counter-jet channel of the Crab PWN (jBietenholz et al.l 12004 ). Com- 
plex interactions between the PWN and the surrounding ejecta filaments are als o seen in such 



a reg i on, roughly corresp onding to the inner portion of the Norteastern "bay" (IMichel et al. 



199ll ; iHester et al.lll995l ). We retrieved and inspected HST/WFPC2 images of the Crab 
collected in 1994 and in 2001 through the F54 7M filter. We found the outer feature to show 
a variability consistent to the one reported by lHester et al.l (119951 ). corresponding to a local 
surface brightness variation of order 25%, very similar to the hot spot of PSR B0540-69. 
We also note that, rescaled at the LMC distance, the variability in the inner nebula (wisps 
and torus) would be difficult to detect, while the variability of the outer structure would be 
outstanding. 

Coming back to PSR B0540-69, the hot spot lies in the region of the PWN tentatively 
identified as an equatorial torus by iGotthelf fc Wang! (120001 ) . The large variability of the 
feature, coupled to the PWN asymmetry with respect to the pulsar position, as well as 
the comparison with the case of the Crab, may point to an alternative scenario in which 
the northeast-southwest axis corresponds to the direction of a pulsar jet/counterjet. This 
may also be supported by the observation that other PWNe do show (in X-rays) the largest 
variability along the jet direction, with apparent complex motion of bright blobs, although 



on sh orter time scales (e.g. for PSR B1509-58 and Vela. iDelanev et al 



2006; Pavlov et al. 



20031 ) as well as on a smaller physical scale (in Vela, iPavlov et al.l 120031 ) . However, no firm 
conclusions may be drawn based on current data. 



3. Conclusions 

With the discovery of significant variations in its PWN emission, PSR B0540-69 shares 
one more characteristic with the Crab pulsar. Moreover, our multi-epoch study of PSR 
B0540-69 yielded a new assessment of the pulsar proper motion, setting an upper limit of 
290 km s _1 . 

Upcoming WFPC2 observations of PSR B0540-69 will allow us to monitor the PWN 
morphology, while lowering the measurable velocity to ~ 220 km s" 1 . The new HST data 
will be very important to shed light on the geometry and dynamic of the system. High- 
resolution HST polarimetric mode observations, to be collected in the same program, will 
offer unvaluable clues in order to understand the overall structure of the PWN and its 
complex interaction with the cage of filamentary ejecta. PSR B0540-69 will possibly become 
a unique extra-galactic laboratory to study and understand the variability and evolution of 
young PWN systems. 
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Fig. 1.— The field of PSR B0540-69 as observed with HST/WFPC2 through the F555W 
filter in 1995 (600 s) and in 2005 (480 s). Pixel size is 0'.'0455. The pulsar optical counterpart, 
as well as the brightest feature of the PWN (the "Hot Spot") are marked. A grid (20 pixel 
spacing, corresponding to 0'.'91) is overplotted to better visualize the apparent displacement 
of the "Hot Spot". 
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Fig. 2.— The inner field of PSR B0540-69 as observed in 1995 through the F555W filter 
is shown, aft er adaptive gaussian smoothing of the HST/WFPC2 image. The asmooth 
algorithm by lEbeling et al.l (120061 ) has been used, setting the minimum S/N threshold to 
10 and using a maximum smoothing kernel of 10 pixels. Isophotal contours to mark the 
hot spot peak (corresponding to 99%, 95% and 99% of the maximum of the PWN surface 
brightness) are plotted in white. Contours generated using the same criteria for the 1999 
and 2005 images are overplotted. The displacement of the hot spot is apparent. 



